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ABSTRACT 

We extend our calculation of physical parameters of GRB jets by modeling the broadband emission of the 
afterglows 970508, 980519, 991208, 000926, 000418, and 010222. Together with 990123, 990510, 991216, 
and 000301c, there are ten well-observed afterglows for which the initial opening angle of the GRB jet can be 
constrained. The jet energies (after the GRB phase) obtained for this set of afterglows are within one decade 
around 5 x 10 50 erg. With the exception of 000418, which requires a jet wider than 1/2 radians, the jet initial 
half-angle in the other cases ranges from 2° to 20°. We find that, in half of the cases, a homogeneous ambient 
medium accommodates the afterglow emission better than the wind-like r~ 2 profile medium expected around 
massive stars. The two types of media give fits of comparable quality in four cases, a wind medium providing a 
better description only for 970508. The circumburst densities we obtain are in the 0.1 — 100 cm" 3 range, with the 
exception of 990123, for which it is below 10~ 2 cm" 3 . If in all ten cases the observed GRB durations are a good 
measure of the ejecta deceleration timescale, then the parameters obtained here lead to jet initial bulk Lorentz 
factors between 70 and 300, anticorrelated with the jet initial aperture, and jet masses ar ound 10 _6 M Q . Our 
results on the jet energy, opening Lorentz factor, and evacuation of material until break-out provide constraints on 
theoretical models of GRB jets. 

Subject headings: gamma-rays: bursts - ISM: jets and outflows - methods: numerical - radiation mechanisms: 
non-thermal - shock waves 



1. INTRODUCTION 

The localization of Gamma-Ray Bursts (GRBs) to within a 
few arc-minutes by the Italian-Dutch satellite BeppoSAX, the 
Interplanetary Network, and the Rossi-X-ray Transient Ex- 
plorer have enabled us to carry out ground-based follow-up 
searches for afterglow emission. The current database of multi- 
wavelength (radio, millimeter, optical, and X-ray) observations 
allows us to begin a statistical study of the physical properties 
of GRB afterglows. 

This is a third in a series of papers modelling the broad- 
band emission of GRB afterglows, with the aim of determin- 
ing the total energy in the relativistic ejecta, the jet opening 
angle, the density and profile of the medium in the immedi- 
ate vicinity (^ 10 18 cm) of the burst, and the microphysi- 
cal shock parameters. The first paper (PK01) has presented 
the modelling of the afterglows 990123, 990510, and 991216 
while the second (Panaitescu 2001) analyzed the peculiar after- 
glow 000301c, whose emission fall-off exhibited a sharp break. 
Here we present our results for the afterglows 970508, 9805 19, 
991208, 000926, 000418, and 010222. With the exception of 
000418, after 1 day the decay of their optical emission is steep 
or exhibits a steepening, as expected if the GRB ejecta are well 
collimated (Rhoads 1999), thus their modelling allows the de- 
termination of the jet aperture and energy. Section §|| outlines 
the model used to fit the broadband emission of these after- 
glows. The jetproperties inferred for individual afterglows are 
presented in §B] and the results for the entire set are analyzed in 



2. THE AFTERGLOW MODEL 
2.1. Model Features 

The calculation of the afterglow emission is carried out in the 
standard framework of relativistic ejecta decelerated by an ex- 
ternal medium (Meszaros & Rees 1997), with allowance for the 
effects due to collimation (Rhoads 1999). The equations gov- 
erning the dynamics of jet-medium interaction and those for 
the calculation of the synchrotron and inverse Compton emis- 
sion are presented in KP00, PK00 and PK01. Similar analytical 
treatments of jet dynamics and/or emission of radiation can be 
found in Waxman (1997), Granot, Piran & Sari (1999), Gruzi- 
nov & Waxman (1999), Wijers & Galama (1999), Chevalier & 
Li (2000), and Sari & Esin (2001). The effect of interstellar 
scintillation on the radio afterglow emission (Goodman 1997) 
is taken into account following the treatment of Walker (1998). 

In our treatment, the afterglow modelling has the following 
features: 

1) the jet is considered uniform, with an energy per solid angle 
independent of direction within the jet; 

2) the shocked gas internal energy density is assumed uniform; 

3) the jet dynamics is calculated by following the evolution of 
its energy (which decreases due to radiative losses), mass (in- 
creasing, as the jet sweeps-up the surrounding medium), and 
aperture (which increases due to jet expansion in the comoving 
frame). The equations employed are accurate in any relativistic 
regime; 

4) the afterglow emission is calculated by integrating over the 
jet dynamics the synchrotron and inverse Compton radiation, 
taking into account the spread in the arrival time of photons 
emitted at a given radius. The Compton parameter and its evo- 
lution are calculated from the electron distribution; 
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5) the shock-accelerated electron distribution is assumed a 
power-law 7V(7) oc 7~ p in the random electron Lorentz fac- 
tor 7, starting from a minimum ji up to a high energy break 7*; 

6) the fractional energy in electrons and magnetic field are con- 
stant throughout the jet deceleration; 

7) in general, the observer is assumed to lie on the jet axis. For 
our jet model, observer off-sets lower than the jet initial aper- 
ture produce insignificant changes in the resulting light-curve 
(see Granot et al. 2002), and therefore in the fitting parameters. 

2.2. Model Parameters and Their Determination 

The model has three parameters that give the jet dynamics: 
the initial jet energy Eq, initial half-angle 6q, and external par- 
ticle density n (or the constant A for a wind-like density pro- 
file^] n(r) = Ar~ 2 ), and three parameters related to the mi- 
crophysics of shocks: the fraction eb of the post-shock energy 
density in magnetic fields, the fractional energy e e in electrons 
if they all had the same Lorentz factor 7,, and the power-law 
index p. For p < 2, there are two additional parameters: the 
fractional energy e of the electrons between 7,; and 7,, which 
parameterizes 7*, and the power-law index q > 2 of the elec- 
tron distribution above 7* (for simplicity the cut-off above 7, is 
assumed to be a power-law). In some cases, the curvature of the 
optical spectrum or consistency between the optical and X-ray 
afterglow emission indicates a significant dust extinction in the 
host galaxy. This is taken into account by assuming an SMC- 
like reddening curve and adding an extra model parameter, the 
Ay extinction in the host rest-frame. We note that the initial jet 
Lorentz factor has little effect on the afterglow emission and is 
not considered a free model parameter. 

The spectrum of the afterglow synchrotron emission (Sari, 
Piran & Narayan 1998) has breaks at the self-absorption fre- 
quency v a , injection frequency Vi corresponding to the min- 
imum electron 7$, cooling frequency v c corresponding to the 
electron Lorentz factor for which the radiative timescale equals 
the dynamical time, and cut-off frequency v* associated with 
7*. Generically, the afterglow emission F u can be written as 



F v = F, 



P "a 



(1) 



where the exponent (3 is non-zero for all break frequencies 
between the observing frequency v and the frequency v p — 
min(i/ c , Vi) at which the synchrotron spectrum peaks, F p be- 
ing the flux at v p . For a relativistic spherical outflow and an 
adiabatic expansion, the break frequencies and peak flux are 
given by: 
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As illustrated in equation (jl|), the afterglow light-curve at a 
given frequency is given by the evolution of F p , v a , i/j, v c , and 
v* which, for constant parameters eb, £e and e, are determined 
by the Lorentz factor Y of the jet, its radius r, and external den- 
sity profile n(r). For a relativistic jet and negligible radiative 
losses, conservation of total jet energy leads to Y oc <~ 3 / 8 and 
r oc t x l A in the case of a homogeneous medium and Y oc t > , 
r oc t 1 / 2 for a wind external medium, before the time 



tj = 0.4 (z + 1) (£o,5ono 1 ) 1/3 d ay 



(3) 



when the jet transits between a quasi-spherical expansion and 
a lateral spreading dominated one |J In equation (Q) i?o,50 the 
initial jet energy measured in 10 50 erg and #o,-i the initial jet 
half-opening measured in 0.1 radians. After tj the jet dynamics 
is described by Y oc t~ x l 2 and (to "zeroth order") r ~ constant. 
The resulting evolution of the afterglow spectral characteristics 
in these two asymptotic regimes are summarized in Table 1, to- 
gether with the afterglow light-curve t~ a at frequencies above 
for slowly cooling electrons, i.e. V{ < v c , and for Comp- 
ton parameter Y below or above 1, i.e. electron cooling due to 
synchrotron losses or up-scatterings, respectively. 

Because the afterglow light-curve decay depends only on the 
index p of the electron distribution (or q above 7*), this index 
can be easily determined from observations one knows the 
location of Vi and v c relative to the observing frequency. For 
measurements made more than a few hours after the GRB, the 
injection frequency is below the optical domain, thus the only 
uncertainties are related to v c . Consistency between the decay 
indices a(p) given in Table 1 and the slope (3(p) of the syn- 
chrotron power-law optical spectrum F v oc v~P, where 



1) for Vi < v < v c 



for v c < v , (4) 



(2) 

where z is the afterglow redshift, £0.53 is the fireball energy in 
10 53 erg, no the external medium density y in cm~ 3 , E e and 
Eb have been normalized to 0.1 and 0.01 respectively, td is the 
observer time measured in days, and £>l,28 is the burst lumi- 
nosity distance measured in 10 28 cm. Similar equations can be 
derived for spreading jets, non-relativistic GRB remnants, or 
non-adiabatic expansion. 

1 This constant is proportional to the ratio between the mass loss rate of the star which ejects the wind and the speed of this wind. We denote by A* the value of 
the constant A relative to that corresponding to 10 — 5 Mq ejected per year at a speed of 1000 km/s. 

2 This is also the time when the jet Lorentz factor equals the reciprocal of the jet half-angle, thus the radiation emitted from the jet edge is no longer relativistic 
beamed away from the observer. 



is commonly used to determine both p and the location of v c 
relative to the optical domain. 

Finding the remaining parameters, five if the high frequency 
cut-off v* is above the highest observing frequency, seven in the 
opposite case, is conditioned by the localization of the spectral 
breaks at some time (not necessarily the same for all breaks), 
either from the afterglow flux at two frequencies bracketing a 
given break, or from the passage of that break through an ob- 
serving band. If we know v a , vj. v c and F p from observations, 
inverting the set of equations (g) above allows the calculation 
of £q, n, £ e and eb- Then, if the afterglow decay exhibited a 
steepening which can be identified as the "jet-break", equation 
([3]) gives the initial jet aperture #o- Together with £q, this allows 
the calculation of the jet initial energy Eq. 

In general the above method cannot be readily used to other 
afterglows because the locations of v a and v c are not suffi- 
ciently constrained by the available data. For instance, evidence 
for self-absorption at radio frequencies exists only for the after- 
glows 970508, 991208, 000301c and, perhaps, 991216. Fur- 
thermore, the approximations usually made in analytical treat- 
ments of the afterglow emission (e.g. Waxman 1997, Wijers 
& Galama 1999, PK00, Sari & Esin 2001) are accurate only 
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after the jet-break time, suggests that the tj 22 100 days, if the 
assumption of relativistic motion is still correct. 

In order to have an observational constrain on the jet initial 
opening, we make the assumption that the brightening seen at 
1 day is due to an observer location 8 {, s which, initially, is 
outside the jet (i.e. o b s > 9o). In this case, the observer will 
see a rising light-curve when the jet has decelerated down to a 
Lorentz factor T ~ (0 o bs — for which the radiation from 
the nearest part of the jet becomes visible to the observer. We 
also assume that the GRB and the dimmer afterglow emission 
prior to the brightening were due to some less energetic ejecta 
located outside the main jet and moving toward the observer, 
but ignore their effect on the dynamics of the jet and afterglow 
emission after 1 day. 

The parameters of the best fit to the radio (Frail et al. 2000b), 
millimeter (Bremer et al. 1998), optical (Sahu et al. 1997, 
Galama et al. 1998a, Pedersen et al. 1998, Sokolov et al. 1998, 
Zharikov et al. 1998), and X-ray (Piro et al. 1998) emission of 
970508 after the onset of the brightening phase, obtained with 
a homogeneous external medium, are given in Table 2. The 
observer is located at o b s ~ 4/3 6q. The fit is not satisfac- 
tory (thus we do not determine confidence intervals), exhibit- 
ing a slower brightening and weaker X-ray emission than ob- 
served. The jet isotropic equivalent energy £q derived by Wijers 
& Galama (1999) for the 970508 afterglow from its spectrum 
at 12 days is twice smaller than our value^ while their value for 
eb is twice larger. The largest discrepancies are in n, our value 
being 25 times larger than theirs, and in e e for which we find 
a value five times larger. From the same 12 day spectrum of 
970508, Granot et al. (1999) found £ 15 times smaller, n four 
times larger, and e b three times smaller than our values, and a 
similar e e . Through an analytical treatment of the radio emis- 
sion of 970508, Frail et al. (2000b) inferred a jet energy 4 times 
smaller than obtained by us, electron and magnetic parameters 
close to equipartition and n ~ 1 cm -3 , close to our results. 

A significantly better fit can be obtained for a wind-like ex- 
ternal medium. This fit, shown in Figure 1, has Eq = 1.6 x 10 51 
erg, O = 18°, 9 obs ~ 5/3 9 , A* = 0.39, e e = 0.15, 
e B = 0.10, p = 2.32, and \ 2 = 570 for 279 degrees of free- 
dom (dof). Although it describes well the rise of the optical 
emission, it does not accommodate the X-ray emission dur- 
ing the brightening. Using a spherical model, Chevalier & Li 
(2000) estimated similar values for A* and Eb, and £q and n 
twenty and seven times smaller, respectively. 

The spectral properties of the best fit models discussed above 
are similar. In both cases the injection frequency vi is slightly 
below the optical domain at 1 day^ , its evolution yielding the 
spectral softening observed by Galama et al. (1998b). At < 100 
days, Vi passes through the radio domain, as inferred from ob- 
servations by Frail et al. (2000b). 

3.2. GRB 980519 

The optical emission of this afterglow had a break of mag- 
nitude A a ~ 0.5 at t ~ 1 day, with a temporal index a Q = 
2.22 ± 0.04 (Jaunsen et al. 2001) after the break, close to that 
measured in X-rays, a x = 2.25 ± 0.04, at about 1 day (Nicas- 
tro et al. 1999). The equality of these two indices is consistent 
with the achromatic break expected in the jet model. At < < 1 
day the slope of the optical spectrum dereddened for Galactic 
extinction, (3 Q = 1.20 ± 0.25 (Halpern et al. 1999), is shallower 

3 For the e e and eg given for 970508 in Table 2, radiative losses are significant (80% until around 10 days), thus the jet energy we infer is likely to exceed that 
obtained by other researchers using adiabatic models for the jet dynamics 

4 Due to the off-axis location of the observer, the passage of Vi through the optical occurs later than predicted by eq. (ph 



over a limited time interval. Various departures from those 
approximations, such as 1) moderately relativistic jets, with 
r of several, 2) jets transiting between quasi-collimated and 
lateral-spreading expansion, 3) electron radiative cooling not 
dominated by a single emission process (synchrotron or inverse 
Compton), 4) afterglow spectral breaks smoothed by the dif- 
ferential relativistic boost and photon arrival-time over the jet 
surface, 5) time changing ordering of the spectral, during the 
afterglow evolution, require numerical calculations to yield a 
more reliable determination of jet parameters. 

3. COLLIMATED AFTERGLOWS 

The model outlined above was used to model the broadband 
emission of ten afterglows to determine the parameters Eq, 6q, 
n (or A* for a wind), e e , eb, and p (plus Ay, e and q in those 
cases where they are relevant) by x 2 -minimization. Nine of 
these afterglows (970508, 980519, 990123, 990510, 991208, 
991216, 000301c, 000926, and 010222) were selected based 
on the existence of a break in (or a steep decay of) the optical 
light-curve, allowing the calculation of the jet initial opening, 
and sufficient broadband observations to make the modelling 
meaningful. The 000418 afterglow has been added due to its 
good multiwavelength coverage, although its emission does not 
exhibit a signature of a jet-break. 

In calculating the afterglow optical fluxes we assumed a 5% 
error in the magnitude-to-flux conversion and Galactic redden- 
ing and we subtracted the reported contributions of the host 
or contaminating galaxies. X-ray fluxes have been calculated 
from the reported band fluxes (2-10 keV, usually) and X-ray 
spectral slopes. 

The modelling of the afterglows 990123, 990510, 991216, 
and 000301c is presented elsewhere (PK01, Panaitescu 2001). 
The use of a larger set of X-ray measurements for the 990123 
afterglow (Costa 1999) and the reduction of the assumed er- 
ror in the Galactic extinction of the optical emission of 991216 
from 10% to 5%, led in these two cases to small changes in 
the best fit parameters. Below we discuss the features of the 
broadband emission and the modelling of the remaining six af- 
terglows. 

3.1. GRB 970508 

During the first day after the burst, the optical emission of 
970508 exhibits a brightening by more than one magnitude, 
lasting for about 1 day, followed by a long-lived decay of tem- 
poral index a a = 1.17 ± 0.03 (Sokolov et al. 1998). Simul- 
taneous with the brightening, Galama et al. (1998b) observe 
a reddening of the optical spectrum, the slope changing from 
f3 ~ 0.5 at 1 day to f3 ~ 1.1 at t > 2 days. They interpret the 
softening of the optical spectrum as the passage of the cooling 
frequency i/ c , thus the latter (3 Q implies p = 2.2 . 

The radio emission of this afterglow has been monitored for 
about 1 year. Frail, Waxman & Kulkarni (2000b) report that the 
slope of the radio spectrum evolved from [3 r = — 0.25 ± 0.04 
before 80 days to an average j3 r = 0.6 ± 0.2 after 100 days, 
which must be due to the passage of the injection frequency 
Therefore the latter r implies that p = 2.2 ± 0.4, as the 
radio domain is expected to lie below v c . After 100 days, the 
radio emission decayed as a power-law of index a r = 1 .2 ± 0.1 
(Frail et al. 2000b), which, given the above p and that a ~ p 
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than that measured by Nicastro et al. (1999) at about the same 
time in X-rays, (3 X = 1.72 ± 0.42. The difference between the 
two slopes is close to that expected when v c is between opti- 
cal and X-rays but, given the their large uncertainties, does not 
provide a compelling proof. 

Numerically we find that the radio (Frail et al. 2000a), opti- 
cal (Vrba et al. 1999, Jaunsen et al. 2001) and X-ray (Nicastro 
et al. 1999) emission of 980519 can be well accommodated by 
a spreading jet interacting with a homogeneous medium, and 
with v c between optical and X-rays (Figure 2). The best fit 
with a jet interacting with wind medium has Eq = 1.1 x 10 51 
erg, O = 6.7°, A* = 3.5, e e = 0.036, e B = 0.22, p = 2.43, 
and x 2 — 73 for 46 dof. This model yields a shallower break 
than observed in the /-band light-curve of 9805 19 and provides 
a poorer fit to the radio data. 

3.3. GRB 991208 

The radio emission of this afterglow (Galama et al. 2000) 
exhibited a quasi-flat behavior until ~ 10 days (Figure 3), fol- 
lowed by a power-law decay (Galama et al. 2002) which is 
much shallower than the £~ 2 - 2±0 - 2 observed in the optical at 
2-7 days (Castro-Tirado et al. 2001). The flatness of the early 
radio emission indicates that the V{ frequency is above the ra- 
dio domain and that either the external medium is homogeneous 
and the GRB remnant is a jet observed after tj (i.e. after the jet- 
break) or the ambient medium is wind-like and the GRB jet is 
seen before tj. This suggests that the steepening of the radio 
emission at 10 days is due to the ^-passage. Then the shallow 
radio decay after 10 days requires a hard electron distribution, 
p < 2, and the steep optical decay implies that the -break is 
below the optical domain. 

The best fit to the data obtained with a jet interacting with 
a homogeneous medium is shown in Figure 3 and has the pa- 
rameters given in Table 2. The spectral characteristics (break 
frequencies and peak flux) of the model afterglow emission are 
consistent with those obtained by Galama et al. (2000) by fit- 
ting the spectrum of 991208 at four epochs. The electron dis- 
tribution cut-off at higher energy is characterized e = 0.49 and 
q = 2.7. 

The data can be fit equally well (x 2 = 110 for 97 dof) 
with wind medium of A* = 0.65 and a jet with parameters 
E = 3.2 x 10 50 erg, 6 = 14°, e e = 0.054, e B = 0.021, 
p ~ 1.4, e = 0.32, and q — 3.1. The last five parameters 
above and the implied £q are similar to those determined by Li 
& Chevalier (2001). 

3.4. GRB 000418 

The i?-band emission of this afterglow exhibited a flattening 
after only 10 days (Klose et al. 2000, Berger et al. 2001), in- 
dicating an underlying host galaxy. By fitting the optical light- 
curve with the power-law fall-off expected of GRB afterglows 
plus the constant contribution of the host, Berger et al. (2001) 
obtain a decay index a a = 1.41 ± 0.08. The 8.5 GHz radio 
emission exhibits a gradual steepening to a power-law of simi- 
lar index, a r = 1.37±0.10(Bergeretal. 2001) after - 40 days, 
which, most naturally, is due to the v\ frequency falling below 
the radio domain. The emission of the 000418 afterglow does 
not exhibit a steepening that can be attributed to the spreading 
of a jet, therefore its modeling should not constrain significantly 
the initial jet aperture. 

Optical observations at other frequencies are rather scarce, 
making the determination of the optical spectral slope rather 



difficult. From the only two simultaneous K and R measure- 
ments at times when the afterglow emission is dominant, we 
infer f3 Q = 1.62 ±0.15 (note that Klose etal. 2000 find a redder 
spectrum, with [3 = 1.90 ± 0.15). Then consistency between 
a (Table 1) and (3 (eq. [[jj) suggests that the cooling frequency 
is below the optical range and/or there is a significant dust red- 
dening within the host galaxy. 

With our jet model, the best fit for this afterglow is obtained 
with a wide jet (8q ~ 1 radian) interacting with either a homo- 
geneous medium or a wind. Table 2 lists the parameters in the 
former case, the model light-curves being shown in Figure 4. In 
the latter case, the best fit parameters are Eq — 2.2 x 10 51 erg, 
O = 60°, A* = 0.69, e e = 0.10, e B = 0.027, and p = 2.04. 
The two models have x 2 — 55 and \ 2 — 56, respectively, for 
61 dof. Taking the host extinction Ay as a free parameter, does 
not improve the fits significantly. Compared to the parameters 
inferred by Berger et al. (2001), the above n and A* are larger 
by a factor 100 and 10, respectively, e e is twice larger, e B is 
smaller by a factor 10-30, and the isotropic £q for the Eq and 
6*o above is smaller by a factor of a few. We note that jets nar- 
rower than 1/3 radians yield poorer fits, with x 2 *J 130 for 61 
dof, while acceptable fits can be obtained with spherical ejecta. 

3.5. GRB 000926 

The X-ray emission (Piro et al. 2001) of this afterglow pro- 
vided for the first time evidence (Harrison et al. 2001) that the 
X-ray emission may be inverse Compton scatterings (PK00, 
Sari & Esin 2001). This is suggested by that the extrapolation 
of the optical spectrum, after dereddening for the host (intrin- 
sic) extinction, falls below the observed X-ray flux. 

The optical emission of 000926 exhibited a break of magni- 
tude Aa ~ 0.75 at few days, with a post-break temporal index 
a a ~ 2.35 ± 0.05 (Fynbo et al. 2001, Price et al. 2001). If 
interpreted as a jet break, it requires that p 2.4, which would 
imply (eq. [f|]) an optical spectrum significantly harder than ob- 
served at t ~ 1 day: j3 = 1.42 ± 0.06 (Fynbo et al. 2001) or 
P = 1.53 ± 0.07 (Price et al. 2001). Within the fireball model, 
consistency between the optical spectral slope and temporal in- 
dex requires a significant host extinction. From the curvature 
of the near infrared-optical spectrum, Fynbo et al. (2001) in- 
fer Ay — 0.18 ± 0.06, corresponding to an extinction in the 
observer /-band of 0.4 ±0.1 magnitudes, thus the dereddened 
afterglow spectrum has an optical slope j3 a ~ 1. Then equation 
(Q) and p ~ 2.3 imply that v c is below the optical domain. 

For a homogeneous medium, the best fit obtained with a 
model with the above features (Figure 5) has parameters (table 
2) that are close to those obtained Harrison et al. (2001), except 
Eo and e B , for which we find values three times smaller and 
eight times larger, respectively. The best fit model with a wind 
medium has x 2 — 270 for 102 dof, yielding radio fluxes larger 
than observed, and parameters Eq = 2.7 x 10 51 erg, 6q = 2.0°, 
A* = 2.0, e e = 0.042, e B = 1.6 x 10~ 4 , andp = 2.70. Note, 
however, that Harrison et al. (2001) found a significantly better 
fit (X 2 — 167 for 1 14 data points) for a wind medium. 

3.6. GRB 010222 

The index a of the power-law decay of the optical emission 
of 010222 steepened by Aa = 0.6 ± 0.1 at about 0.5 days, 
to a Q = 1.30 ± 0.05 (Masetti et al. 2001, Stanek et al. 2001), 
a value similar to that seen in the X-rays, a x = 1.33 ± 0.04 
('t Zand et al. 2001). The jet interpretation of this break re- 
quires an electron distribution with p ~ 4/3. As in the case 
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of 991208, a hard electron distribution lowers the cut-off fre- 
quency v* sufficiently to yield a break of the afterglow de- 
cay when v* passes through of the observing band. This pro- 
vides a natural explanation for the second steepening observed 
in the optical after 10 days by Fruchter et al . (2001) (see also 
http://www. stsci.edu^fruchter/GRB/0 1 0222j ). 

The low index p required by the jet interpretation of the first 
break also implies an intrinsic optical spectrum harder than 
p/2 ~ 2/3, i.e. harder than observed: j3 a = 0.89 ± 0.03 (Jha et 
al. 2001, Lee et al. 2001), indicating the existence of a signif- 
icant dust reddening in the host galaxy. The best fit we found 
to the radio (Berger & Frail 2001), optical (Cowsik et al. 2001, 
Masetti et al. 2001, Sagar et al. 2001, Stanek et al. 2001) and 
X-ray ft Zand et al. 2001) data has Ay = 0.21 and a large 
X 2 — 236 for 87 dof. More than half of this \ 2 arises from 
8 optical and X-ray data, suggesting either that some reported 
observations have underestimated uncertainties or that there are 
short timescale fluctuations in the afterglow emission (Cowsik 
et al. 2001). The jet model presented in Figure 6 faces also 
another problem: even if electrons acquire 100% of the post- 
shock fluid energy, the passage of through the X-ray domain 
takes place too early, at about 1 day, leading to a discrepancy 
between the last few X-ray data and the model expectations. 
Given the unsatisfactory fit provided by this model, we do not 
include in Table 2 the uncertainties of the best fit parameters. 

As shown in Figure 6, the millimeter model emission falls 
below the detections reported by Fich et al. (2001) and Kulka- 
rni et al. (2001) at 220 GHz and 350 GHz, and below the upper 
limit on the 95 GHz emission found by Bremer et al. (2001). 
The constancy of the observed fluxes over almost one decade 
in time and the steep spectrum between 220 GHz and 350 GHz 
are hard to accommodate within the jet model, suggesting that 
the millimeter excess seen in 010222 is due to a dusty, star- 
forming host galaxy (Kulkarni et al. 2001). These data were 
not included in our fits. 

The best fit obtained with a wind-like medium has a similar 
large % 2 = 241 for 87 dof, and parameters E = 2.5 x 10 50 
erg, O = 3.4°, A, = 0.18, e e = 1.4 x 10~ 2 , e B = 1.0 x 10~ 3 , 
p = 1.43, and Ay = 0.22. It provides a better description of 
the last X-ray measurements than the homogeneous medium 
model, but a poorer one for the decay steepening seen in the op- 
tical at 0.5 days, as jets interacting with winds yield long-lived, 
smooth light-curve breaks (KP00). The millimeter emission in 
the wind model is also consistent with the 95 GHz upper limit 
and falls below the 220 GHz and 350 GHz detections. 

4. JET PROPERTIES 

4.1. Jet Energy 

For eight of the ten GRBs jets whose basic parameters are 
listed in Table 2, the ejecta kinetic energy Eq at the beginning 
of the afterglow phase is between 10 50 and 5 x 10 50 erg. A 
narrow distribution of the jet kinetic energy has also been in- 
ferred by Piran et al. (2001) for a larger set of afterglows, based 
on the width of the observed X-ray luminosity distribution at 
0.5 days. Note, however, than the jet energy for the 970508 and 
000418 afterglows is significantly larger, being around 2 x 10 51 
ergs. 

The initial half-angle of the jet 0q is correlated with Eq: ex- 
cluding 000418, for which we find Eq and 8q much larger than 
for the other afterglows, their linear correlation coefficient is 

5 No value is given for GRB 970508, as we interpreted its afterglow emission 
that has produced the GRB emission 



r(E , ) = 0.68 ± 0.04, i.e. a 5% probability of obtaining this 
correlation by chance. 

The 7-ray energy output E 1 for our set of afterglows, calcu- 
lated from the fc-corrected isotropic-equivalent £ 1 of the GRB 
emission in the 20-2000 keV obtained by Bloom, Frail & Sari 
(2001) and the jet aperture resulting from afterglow modeling, 
i.e. E 7 = (1/2)£ 7 (1 - cosflo), are listed in Table 3. We find 
that, excluding 000418, for which the jet opening is much larger 
than for the other cases analyzed here, the dynamical range 
of E 1 is 30, i.e. a factor 3 larger than obtained by Frail et 
al. (2001) after determining the jet aperture from the jet-break 
time (eq. 

Table 3 also shows the resulting 20-2000 keV GRB effi- 
ciency^] defined by e 7 = E 1 j(E 1 + E ). Note that, with the 
exception of 980519, e 7 is larger than ~ 50%, most likely ex- 
ceeding the ability of internal shocks to channel the dissipated 
energy into the 20-2000 keV band. This suggests that, dur- 
ing the GRB phase, jets are inhomogeneous on angular scales 
smaller than 8q, such that we are biased toward observing bursts 
whose outflow has a bright spot moving directly toward the ob- 
server or very close to the observer's line of sight toward the 
GRB source (Kumar & Piran 2000). Then the true efficiency 
of the GRB is much smaller and the jet energy Eq at the begin- 
ning of the afterglow phase is closer to the jet energy before the 
GRB phase (Piran et al . 200 1 ) . 

4.2. External Medium 

Our results show that models with a homogeneous medium 
can accommodate the broadband emission of all EIGHT af- 
terglows, while wind-like medium is consistent with the ob- 
servations in at most four cases. If our assumption regard- 
ing the jet uniformity is correct, then a GRB model involving 
a massive star is allowed in the remaining four cases only if 
there is a mechanism for homogenizing the wind surrounding 
the star prior to its interaction with the jet. Ramirez-Ruiz et 
al. (2001) have shown that the interaction between the wind of 
a Wolf-Rayet star and a circumstellar medium of n = 1 cm -3 
leads to the formation of a quasi-uniform, hot shell of density 
~ 10 3 cm~ 3 , extending from > 10 16 cm up to ~ 10 18 cm. 
More tenuous (or colder) media could produce thicker and less 
dense shells, consistent with the range of densities found here. 

The particle density given in Table 2 for homogeneous media 
range from values typical for the interstellar medium (970508, 
980519, 990510, 010222) to those of diffuse hydrogen clouds 
(991208, 991216, 000301c, 000418, 000926). In one case 
(990123) we find an external density below 10~ 2 cm~ 3 , char- 
acteristic of a hot component of the interstellar medium or a 
galactic halo. A similar low density was also obtained for the 
afterglow 980703 (PK01). These values are 2-5 orders of mag- 
nitude smaller than those implied by the Njj column densi- 
ties inferred by Galama & Wijers (2001) for 970508, 980519, 
980703, 990123, and 990510 from the absorption seen in their 
soft X-ray spectra. Furthermore, external densities higher than 
inferred by us are expected if GRBs are related with the death 
of massive stars, as in the collapsar model (Woosley 1993, 
Paczynski 1998, MacFadyen & Woosley 1999). Our results are 
compatible with the above results/expectations if the gas in the 
vicinity of the GRB was evacuated prior to the jet ejection. Re- 
cently Scalo & Wheeler (2001) have pointed out that the super- 
novae and H II region winds occurring in a cluster of massive 
stars form "superbubbles" within giant molecular clouds, with 

ng from a jet seen from outside its initial opening, i.e. not from the same ejecta 
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local densities that range over few orders of magnitude, possi- 
bly being as low as 1CP 3 cm~ 3 , depending on the superbubble 
age, ambient medium and power input from supernovae. 

4.3. Jet Lorentz Factor and Mass 

The afterglow emission is only weakly dependent on the ini- 
tial jet Lorentz factor To, which determines the evolution of the 
radiative losses in the early afterglow. Thus Tq cannot be di- 
rectly constrained through afterglow modelling. However, the 
inferred jet parameters can be used to determine the jet Lorentz 
factor T during the afterglow phase: 



400 



0,50 



9 0,-l«0 



1 + z 



-3/8 



(5) 



where the usual notation X — 10 n X n was used and t is mea- 
sured in seconds. Thus Tq can be calculated if one knows when 
the afterglow began, i.e. the jet deceleration timescale to- In a 
few bursts (Giblin et al. 1999, Tkachenko et al. 2000), soft X- 
ray emission has been observed from the end of the GRB phase 
up to 10 4 s, indicating that the external shock had already set 
in by the end of the GRB. In other cases (Pian et al. 2001, 't 
Zand et al. 2001), no X-ray emission has been detected after 
the GRB, suggesting that to is larger then the burst duration. To 
constrain Tq, we assume^ that the observed GRB duration is a 
good measure of to. Equation (Q) shows that T has a moderate 
dependence on t, thus the error due to this assumption is, likely, 
not too large. 

Table 3 lists the values of To obtained for the best fit parame- 
ters given Table 2. As shown in Figure 7, Tq varies between 70 
and 300 and is anticorrelated with the jet initial opening angle. 
Their linear correlation coefficient r (To, 6*o ) = —0.47 ± 0.11 
corresponds to a 20% chance of obtaining by chance this corre- 
lation in the null hypothesis. The anticorrelation of Tq with 8o 
has been recently suggested by Salmonson & Galama (2002) 
based on the positive correlation they observed in several cases 
between the GRB pulse lag-time and the afterglow jet-break 
time. We note, however, that the dependence they infer between 

To and 9q (To oc Oq 8 ^ 3 ) is much stronger than that found by us 
from afterglow modeling (To oc 0q ' 3 , see Figure 7). 

That wider jets have lower bulk Lorentz factors may also 
be the origin of the GRB pulse lag-time anticorrelation with 
the burst peak luminosity found by Norris, Marani & Bonnell 
(2000) and Salmonson (2000). For the same kinetic energy, nar- 
rower jets have a larger energy per solid angle, which could lead 
to a higher GRB peak luminosity. If, in the comoving frame of 
the burst, the pulse duration and/or peak time dependence on 
photon energy are set by a process (e.g. electron cooling) whose 
timescale is the same in all bursts then, due to the relativistic 
contraction of time, GRBs from faster jets would have smaller 
pulse lag-times. Then the To — Oo anticorrelation implies that 
(narrower) jets with higher peak luminosities (are faster and) 
yield shorter 7-ray pulse lags. 

As indicated in Table 3, the product To^o ranges from 10 
to 80, implying that, during the GRB phase, due to relativistic 
beaming, the observer receives emission from only a small frac- 
tion, less than 2%, of the jet surface. Thus calculations of the 
jet 7-ray output E 1 obtained by equating the energy per solid 
angle in the Tq 1 region visible during GRB phase to the energy 

6 This assumption is also used for GRB 970508 although, in our interpretation 
different ejecta. The inclusion of 970508 does not change the following conclusions 



per solid angle within the much wider region seen during the 
afterglow phase could lead to rather unreliable results. 

From the jet energy Eq at the beginning of the afterglow 
phase, one can also calculate the initial et mass: Mq = 
c~ 2 Eo/To. The results are given in Table 3 and Figure 7. The 
jet mass is correlated with the jet opening, the linear correlation 
coefficient being r{9 , M ) = 0.96 ± 0.03. We note that M jet 
increases slower than 9q, thus the ratio between the jet mass and 
the stellar mass within the jet opening decreases with increas- 
ing 8q. For a 10 Mq GRB progenitor and the jet masses given 
in Table 3, this ratio is between 10~ 5 and 10~ 4 , indicating that 
prior to the jet release the stellar material along the jet direction 
is strongly evacuated. 



4.4. Microphysical Parameters 

The results of Table 2 show that the fractional energy in the 
magnetic field spans three orders of magnitude and that the in- 
dex p of the power-law distribution of shock-accelerated elec- 
trons is not universal. In four of the afterglows analyzed here, 
the shallow fall-off of either the radio or the optical light-curve 
after the jet break requires p ~ 1.5. Meszaros, Rees & Wijers 
(1998) have shown that, for a fixed p, variations in the jet en- 
ergy per solid angle could lead to range of light-curve decay. 
Because the observer receives radiation from the entire jet sur- 
face after the jet-break time tj, the internal structure of the jet 
has little effect on the light-curve decay index after tj, thus we 
believe that the values of p determined by modelling the post 
jet-break afterglow decay are not sensitive to the angular struc- 
ture of the outflow. 

We note that, for a fractional energy electron e close to 
equipartition, the hard electron distributions [p < 2) identified 
in the 991208, 991216, 000301c, and 010222 afterglows, lead 
to a j/* -break passing through the optical band at/after few days, 
yielding the steepening seen in the optical emission of these af- 
terglows. 

5. CONCLUSIONS 

Our modelling of the broadband emission of ten afterglows 
reveals several properties of GRB jets, which represent con- 
straints on the models for GRB progenitors (Woosley 1993, 
Paczyriski 1998, Vietri & Stella 1998, MacFadyen & Woosley 
1999, Meszaros, Rees & Wijers 1999, MacFadyen, Woosley & 
Heger2001): 

1) the jet energy has a relatively narrow distribution, the values 
determined here being within a factor of 5, around ~ 5 x 10 50 
erg, 

2) the jet initial Lorentz factor is between ~ 100 and 300, 

3) narrower jets are less massive and more relativistic than 
wider jets 

4) the baryonic mass encountered by the jet (as it breaks out) 
is less than 10~ 4 of the material that the GRB progenitor had 
initially within the jet aperture, 

5) the surrounding medium does not have, in general, the r~ 2 
profile expected for the unperturbed wind of a massive GRB 
progenitor. In most cases we find that the density of the exter- 
nal medium is between 0.1 cm~ 3 and 100 cm~ 3 . 

The conclusions and the jet parameters presented here were 
obtained by modelling the afterglow data within a specific 
framework and under certain assumptions, the most notable be- 
ing the uniformity of the jet and the constancy of the energy 

of the brightening of its afterglow, the burst and afterglow emission arise from 
regarding the jet initial Lorentz factor and mass 
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release parameters (e e , eb)- For simplicity, the observer was 
located on the jet symmetry axis. Until the time when first ob- 
servations are done (few hours to 1 day), the narrow jets con- 
sidered here undergo significant lateral spreading, so that the 
afterglow light-curves seen by an observer located off the jet 
axis (but still within the initial jet opening, to allow the GRB 
to be detected and localized) differ little from those seen by an 
on-axis observer. 

More complex jet models for GRB afterglows, such as that 
of a structured jet proposed by Rossi et al. (2001), or a hy- 
drodynamical treatment of the jet lateral spreading (Granot et 



al. 2001), may yield different jet parameters and constraints on 
GRB progenitors than presented here. We note that the exis- 
tence of a quasi-universal jet energy has also been established 
in a less model-dependent way by Piran et al. (2001), based on 
the narrow width of the afterglow X-ray luminosity at 1/2 day, 
when a good fraction of the entire jet is visible to the observer, 
thus this property should also be present in more sophisticated 
jet models. 

AP acknowledges the supported received from Princeton 
University through the Lyman Spitzer, Jr. fellowship. 



REFERENCES 



Berger, E. & Frail, D. 2001, GCN #968 

Berger, E. et al. 2001, ApJ, 556, 556 

Bloom, J., Frail, D. & Sari, R. 2001, AJ, 121, 2879 

Bremer, M., Castro-Tirado, A. & Moreno, R. 2001, GCN #995 

Bremer, M. et al. 1998, A&A, 332, L13 

Castro-Tirado, A. et al. 2001, A&A, 370, 398 

Chevalier, R. & Li, Z. 2000, ApJ, 536, 195 

Costa, E. 1999, A&AS 138, 425 

Cowsik, R. et al. 2001, BASI, 29,91 

Fich, M. et al. 2001, GCN#971 

Frail, D. et al. 2000a, ApJ, 534, 559 

Frail, D., Waxman, E. & Kulkarni, S. 2000b, ApJ, 537, 191 

Frail, D. et al. 2001, ApJ, 562, L55 

Fruchter, A. et al. 2001, GCN#1()87 

Fynbo, J. et al. 2001, A&A, 373, 796 

Galama, T. et al. 1998a, ApJ, 497, L13 

Galama, T. et al. 1998b, ApJ, 500, L97 

Galama, T. et al. 2000, ApJ, 541, L45 

Galama, T. & Wijers, R. 2001, ApJ, 549, L209 

Galama, T. et al. 2002, in preparation 

Giblin, T. et al. 1999, ApJ, 524, L47 

Goodman, J. 1997, New Astronomy, 2, 449 

Granot, J., Piran, T. & Sari, R. 1999, ApJ, 527, 236 

Granot, J. et al. 2001, in "GRBs in the Afterglow Era", eds. E. Costa, F. Frontera 

& J. Hjorth, Springer- Verlag:Berlin, 312 
Granot, J„ Panaitescu, A., Kumar, P. & Woosley, S. 2002, ApJ, submitted 

( |istro-ph/0201322| ) 
Gruzmov, A. & Waxman, E. 1999, ApJ, 511, 852 
Halpern, J., Kemp, J., Piran, T. & Bershady, M. 1999, ApJ, 517, L105 
Harrison, A. et al. 2001, ApJ, 559, 123 
Jaunsen, A. et al. 2001, ApJ, 546, 127 
Jha, S. et al. 2001, ApJ, 554, L155 
Klose, S. et al. 200, ApJ, 545, 271 
Kulkarni, S. et al. 2001, GCN #996 
Kumar, P. & Panaitescu, A. 2000, ApJ, 541, L9 (KP00) 
Kumar, P. & Piran, T. 2000, ApJ, 535, 152 
Lee, B. et al. 2001, ApJ, 561, 183 
Li, Z. & Chevalier R. 2001, ApJ, 551, 940 
MacFadyen, A. & Woosley, S. 1999, ApJ, 524, 262 
MacFadyen, A., Woosley, S. & Heger, A. 2001, ApJ, 550, 410 



Masetti, N. et al. 2001, A&A, 374, 382 

Metzger, M. et al. 2000, GCN #733 

Meszaros, P. & Rees, M.J. 1997, ApJ, 476, 232 

Meszaros, P., Rees, M.J. & Wijers, R. 1998, ApJ, 499, 301 

Meszaros, P., Rees, M.J. & Wijers, R. 1999, New Astronomy, 4, 303 

Nicastro, L. et al. 1999, A&A, 138, S437 

Norris, J., Marani, G. & Bonnell, J. 2000, ApJ, 534, 248 

Paczyhski, B. 1998, ApJ, 494, L45 

Panaitescu, A. & Kumar, P. 2000, ApJ, 543, 66 (PK00) 

Panaitescu, A. & Kumar, P. 2001, ApJ, 554, 667 (PK01) 

Panaitescu, A. 2001, ApJ, 556, 1002 

Pedersen, H. et al. 1998, ApJ, 496, 311 

Pian, E. et al. 2001, A&A, 372, 456 

Piran, T., Kumar, P., Panaitescu, A. & Piro, L. 2001, ApJ, 560, L167 
Piro, L. et al. 1998, A&A, 331, L41 
Piro, L. et al. 2001, ApJ, 558, 442 
Price, P. et al. 2001, ApJ, 549, L7 

Ramirez-Ruiz, E., Dray, L., Madau, P. & Tout, C. 2001, MNRAS, 327, 829 
Rhoads, J. 1999, ApJ, 525, 737 
L ag 
ph/01 12083] 
Sagar, R. et al. 2001, BASI, 29, 91 
Salmonson, J. 2000, ApJ, 544, LI 15 



Rossi, E., Lajszati, D, & Rees, M.J. 2002, MNRAS, submitted (fr stro- 



Salmonson, J. & Galama, T. 2002, ApJ, submitted ( astro-ph/0 112298 ) 

Sahu, K. et al. 1997, ApJ, 489, L127 

Sari, R., Piran, T, & Narayan, R. 1998, ApJ, 497, L17 

Sari, R. & Esin, A. 2001, ApJ, 548, 787 

Scalo, J. & Wheeler, J.C. 2001, ApJ, 562, 664 

Sokolov, V. et al. 1998, A&A, 334, 117 

Stanek, K. et al. 2001, ApJ, 563, 592 

Tkachenko, A. et al. 2000, A&A, 358, L41 

Vietri, M. & Stella, L. 1998, ApJ, 507, L45 

Vrba, F. et al. 1999, ApJ, 528, 254 

Walker, M. 1998, MNRAS, 294, 307 

Waxman, E. 1997, ApJ, 485, L5 

Wijers, R. & Galama, T. 1999, ApJ, 523, 177 

Woosley, S. 1993, ApJ, 405, 273 

't Zand, J. et al. 2001, ApJ, 559, 710 

Zharikov, S., Sokolov, V., & Baryshev, Y. 1998, A&A, 337, 356 



TABLE 1. Evolution of spectral parameters and indices of power-law light-curves at v > v\. 



n oc t/tj f„ oc Vi oc 


v c oc 

Y < 1 Y > 1 


F p OC 


-dlni^/dlnf 

Vi < V < V c V c < V 

Y < 1 Y > 1 


r° < 1 t° t- 3 / 2 

r -2 < 1 t -3/5 f -3/2 

r°,r- 2 >1 t- 1 / 5 t- 2 


£-1/2 t (3p-8)/(8-2p) 
t l/2 t (3p-4)/(8-2p) 
t Q £(2p-4)/(4-p) 


t° 

£-1/2 

t-P 


IP - ! \p ~ 5 jp-ik 

3 1 3 1 3 p 
4^ 4 4^ 2 AP 8-2p 

p-2 

P P P~ T=p- 
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TABLE 2. Best fit parameters for a homogeneous medium and 90% confidence levels for ten GRB afterglows. 



GRB 


Eq 

(10 50 erg) 


(deg) 


n 

( cm" 3 ) 


£e 
(lO" 2 ) 








P 


X 2 /dof 


n oc r 2 ? 


970508 


20 


18.3 


0.75 


11 


4.5 x 10 


-2 




2.18 


780/279 


yes 


980519" 


4 1+4-8 


9 q+0.2 
z -°-0.2 


14+ - 32 

u - 1 *-0.03 


H-3 


(3-5±2i) 


X 


io- 


-5 


9 70+0. 07 
z -< 8-0.04 


53/46 


no 


990123 


i c+3.3 
i -°-0.4 




(i.9t;-j) lo- 3 


i3t; 


(7-4t 2 5 3 9 ) 


X 


10" 


-4 


9 9Q+0.05 
z - zo -0.03 


55/56 


no 


990510 


1 4+4.9 


o 1+0.1 
,J - i -0.5 


o.29i8:S 




(5-2t 4 4 2 7 ) 


X 


10" 


-3 


1 S q + 0.18 
l.»O_0 .01 


36/69 


no 


991208 


2 4+ 2 - 8 


12.8±^ 


18±| 2 


c c+2.1 
°-°-0.9 


(3.5±g : ?) 


X 


10" 


-2 


1 r.q+0.03 
1 - oo -0.03 


112/97 


yes 


991216 




9 7+0.4 
z - ' -1.0 


4 7+6-8 
^•'-1.8 


1-41S1 


(i-s±S:^) 


X 


10" 


-2 


1 S6+ - 03 

i -°°-0.03 


42/41 


yes 


000301c 


q q+0.3 


iq 7+O.6 


27l* 


6-2±i;S 


(7.2t?'j) 


X 


10" 


-2 


1 - 4d -0.07 


119/96 


maybe 


000418 


OO + 120 


50l 3 ° 


97+250 
Z ' -14 




(6-6^ 8 7 ) 


X 


10" 


3 


2 04+ - 07 

Z - U4 -0.18 


55/61 


yes 


000926 


q O+0.3 


s-i±8:g 


22^5 


iol 2 2 


(6.511:5) 


X 


10" 


-2 


2 40+ - 01 

z - w -0.02 


159/102 


no 


010222 


5.1 


4.6 


1.7 


0.43 


6.7 x 10 


-5 




1.35 


230/87 


yes 



a redshift unknown, z = 1 was assumed. 
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TABLE 3. Burst properties and jet characteristics inferred from the best fit parameters given in Table 2. 



GRB z 4 6) 4 C) 4 d) r< e) T 9 MjQ 







(10 50 erg) 




(s) 






(1O" 6 M ) 


970508 


0.84 






35 


150 


47 


8.2 


980519 


l(<0 


0.52 


0.11 


40 


250 


10 


0.90 


990123 


1.60 


4.9 


0.76 


100 


300 


11 


0.28 


990510 


1.62 


1.3 


0.48 


100 


140 


7.8 


0.53 


991208 


0.71 


18 


0.88 


60 


68 


15 


2.0 


991216 


1.02 


3.0 


0.73 


30 


150 


7.0 


0.43 


000301c 


2.03 


6.6 


0.66 


10 


160 


38 


1.2 


000418 


1.12 


148 


0.82 


30 


90 


78 


20 


000926 


2.07 


15 


0.82 


25 


130 


19 


1.4 


010222 


1.48 


11 


0.67 


120 


110 


9.2 


2.5 



a assumed 

b fc-corrected jet energy output in the 20-2000 keV band (Bloom et al. 2001) 

c efficiency of 7-ray emission 

d observed GRB duration 

e jet Lorentz factor at end of GRB 

1 jet mass 
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FIG. 1 . — Best fit for the afterglow of GRB 970508, obtained with a wind-like medium. Parameters are given in text. A jet seen by an observer located outside the 
jet initial opening was assumed in order to obtain the brightening observed after 0.65 days. The afterglow emission seen before this time was not included in the fit 
and could arise from some ejecta located outside the "central", more energetic jet. The cooling frequency u c is below the optical range. The injection frequency Vi 
passes through the optical domain slightly before 1 day and through the radio domain at ~ 100 days. The jet Lorentz factor T falls below 2 at ~ 80 days. Dotted 
vertical lines indicate the amplitude of the interstellar scintillation, triangles showing 2<r upper limits. The / and V band fluxes have been multiplied, for clarity, 
by the factors indicated. The host galaxy contribution inferred by Zharikov et al. (1998) (I = 24.13 ± 0.28, R = 25.55 ± 0.19, V = 25.80 ± 0.14) has been 
subtracted. 
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FIG. 2. — Best fit for the afterglow of GRB 980519, obtained with a homogeneous medium. The model, whose parameters are given in Table 2, has u c between 
the optical and X-ray domains. The electron cooling is due mostly to inverse Compton scatterings. Optical data has been corrected for Galactic extinction of 
E(B — V) = 0.267 (Jaunsen et al. 2001). For this afterglow a redshift has not been measured. We have assumed z = 1, a value typical for other GRBs. 
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FIG. 3. — Best fit for the afterglow of GRB 991208, obtained with a homogeneous medium. The jet has the parameters given in Table 2 and tj 2 days, thus 
it undergoes significant lateral spreading and its edge is visible at the time of the first observations. T is below 4 at the time of observations. The u c and the cut-off 
frequency u t are below the optical domain. The optical emission after 10 days (shown with open symbols) exceeds the model expectations, suggesting the existence 
of a supernova contribution (Castro-Tirado et al. 2001). The host galaxy contribution (/ = 23.46 ± 0.49, R = 24.27 ± 0.15, V = 24.55 ± 0.16) was subtracted. 
The triangle is a 2a upper limit on the 250 GHz flux. 
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FIG. 4. — Best fit for the afterglow of GRB 000418, obtained with a homogeneous medium, and parameters given in Table 2. The u c is below the optical range, 
the Compton parameter is above unity at all times, and T < 2 after ~ 10 days. The Vi falls below radio frequencies at ~ 20 days, yielding the observed fall-offs at 5 
and 8 GHz. The host contamination of R = 23.9 ± 0.2 (Metzger et al. 2000) has been subtracted from the reported i?-band magnitudes. 
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FIG. 5. — Best fit for the afterglow of GRB 000926, obtained with a homogeneous medium. The parameters of the jet are listed in Table 2. The v c lies below the 
optical domain and the Compton parameter is slightly below unity at times when the optical observations were made. The model X-ray emission is due to inverse 
Compton scatterings, with a significant contribution from synchrotron. Optical measurements have been dereddened for host extinction with an SMC-like reddening 
curve and A v = 0.18 (Fynbo et al. 2001). The contribution of a nearby galaxy, corresponding to / = 24.50 ± 0.11, R = 25.19 ± 0.17, and V = 26.09 ± 0.16 
(Price et al. 2001) has been subtracted from ground-based measurements. 
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FIG. 6. — Best fit for the afterglow of GRB 010222, obtained with a homogeneous medium and parameters given in Table 2, and host extinction of Ay = 0.21, 
for an assumed SMC-like reddening curve. The steepening seen in the optical emission at about 0.5 days is the jet-break, while that at ~ 10 days is due to the passage 
of the u t spectral break. The electron cooling is due mostly to up-scatterings. the v c is slightly below or within the optical domain. Note that the model millimeter 
emission is below the constant fluxes (attributed to the host galaxy) and upper limits reported by Bremer et al. (2001), Fich et al. (2001), and Kulkarni et al. (2001). 
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FIG. 7. — Jet Lorentz factor at the end of the GRB (see j yt.3| ) and initial mass versus the jet initial opening. Open symbols are for the afterglows 970508 and 010222, 
for which the best fits obtained are not satisfactory and parameter uncertainties were not determined. The power-law fits to all the afterglow parameters (including 
970508 and 010222) illustrate that wider jets are less relativistic and more massive. 



